Giant Intrinsic Spin and Orbital Hall Effects in Sr2M04 (M=Ru,Rh,Mo) 
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We investigate the intrinsic spin Hall conductivity (SHC) and the d-orbital Hall conductivity 
(OHC) in metallic d-electron systems, by focusing on the t2g-orbital tight-binding model for Sr2A/04 
(M=Ru,Rh,Mo). The conductivities obtained are one or two orders of magnitude larger than 
predicted values for p-type semiconductors with ~ 5% hole doping. The origin of these giant Hall 
effects is the "effective Aharonov-Bohm phase" that is induced by the d-atomic angular momentum 
in connection with the spin-orbit interaction and the inter-orbital hopping integrals. The huge SHC 
and OHC generated by this mechanism are expected to be ubiquitous in multiorbital transition metal 
complexes, which opens the possibility of realizing spintronics as well as "orbitronics" devices. 

PACS numbers: 72.25.Ba,74.70.Pq,85.75.-d 



Recently, the spin Hall effect (SHE) has been attract- 
ing considerable interest since controlling spin by apply- 
ing electric field is one of the key issues in the field of 
spintronics [1]. The SHE was originally predicted based 
on the assumption of extrinsic electron-scattering [2, 3]. 
Intrinsic SHE in semiconductors is caused by uniform 
spin-orbit interaction (SOI) in systems and is indepen- 
dent of the presence of impurities. It existence was 
predicted by Murakami et aZ. [4] for p-type semiconduc- 
tors, and by Sinova et al.[b\ for a two-dimensional elec- 
tron gas (2DEG) with a Rashba-type SOI. Motivated by 
these predictions, many studies have been performed to 
demonstrate the prominent disorder effects on the spin 
Hall conductivity (SHC) [6-9], the possibility of quantum 
SHE in insulating phase of semiconductors [10-12], and 
the existence of the orbital Hall effect (OHE). According 
to the result of the ab initio calculation [13], SHC and or- 
bital Hall conductivity (OHC) for p-type semiconductors 
increase with hole number nh, and they reach about 100 
and 10 [he~^ ■ il~^cm~^] respectively, for rih = 0.05/cell. 

Very recently, the so-called inverse SHE has been ob- 
served for spin currents injected into conventional metals 
such as Al and Pt from a ferromagnet [14, 15]. The SHC 
observed in Pt is ~ 240 [he~^ ■ il~^cm^^], which is about 
10^ times larger than the experimental SHC in n-type 
semiconductors reported in ref. [16]. This observation 
raises the question of whether a novel mechanism may 
be responsible for the giant SHE in d-electron systems. 
Moreover, a simple theoretical model that is capable of 
describing the giant SHE inherent in d-electron systems 
is highly desirable. 

In this Letter, we propose a new mechanism for the 
giant intrinsic SHC and OHC that originate from the 
d-orbital degrees of freedom, which is absent in semicon- 
ductors. In a multiorbital system, a conduction electron 
acquires an "effective Aharonov-Bohm phase factor" due 
to d-atomic angular momentum in conjunction with the 
SOI and the inter-orbital hoppings. This mechanism is 
responsible for the large SHE. An intuitive explanation 
based on the tight-binding model, which is given in Fig. 
3, appears to successfully predict the characteristics of 



SHE in d-electron systems. We stress that not only the 
SOI play a significant role in the SHE, but so does the 
phase of hopping integrals characteristic of d-electron sys- 
tems. In contrast, the "Dirac monopole mechanism" [4] 
is appropriate when massless Dirac cone dispersion ex- 
ists in close proximity to the Fermi level, as is the case 
in semiconductors. The present mechanism of SHE also 
differs from that in the Rashba-type 2DEG model due to 
the momentum-dependent SOI [5]. We further show that 
the calculated SHC depends strongly on the longitudinal 
resistivity in the high-resistivity region, even though the 
SHE is of intrinsic origin. 

To demonstrate the mechanism of the intrinsic 
SHC and OHC due to d-orbital degrees of free- 
dom, we adopt the square-lattice tight-binding model 
with (dxz,dyz,dxy)-OThitals, which was first proposed 
to describe the realistic electron-structure of Sr2M04 
(M=Ru,Rh,Mo) [17], and atomic LS coupling as the 
SOI. This model is subsequently referred to as the t2g- 
model. We show that the calculated SHC and OHC are 
one or two orders of magnitude greater than theoretical 
values in highly doped p-type semiconductors [13]. The 
present study strongly predicts that giant SHE and OHE 
are ubiquitous in d-electron systems, and suggests a new 
method to analyze them in various metals. 

The Hamiltonian for the i2£rniodel in a square lattice 
is given hy H = Hk + Hso , where Hk is the kinetic term 
and Hso = J2i 2Aij-Sj is the SOI for M 4d-orbitals. The 
band parameters of the i2c,-model are: (i) the intraorbital 
nearest neighbor (NN) hopping integrals ti — t2{= t) 
and t^ for the xz,yz and xy orbitals, respectively, (ii) 
intraorbital second NN hopping tg for xy orbital, and 
(iii) interorbital second NN hopping t' between the xz 
and yz orbitals. It will be shown later that interorbital 
hopping t' is indispensable for SHE and OHE. 

The matrix elements of Hk are then easily given by 
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-2tcoskx, ^2 



-2t cos ky, ^3 
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-2t3{coskx + cos ky) 



4t'^ cos k. 



cos ky 

+ and g 



— At' sm kx sin ky. Here, the first, second, and third rows 
(columns) of Hk correspond to \xz), \yz), and \xy), re- 
spectively, and 1, 2, and 3 denote orbitals xz, yz, and 
xy, respectively. We set h — 1 during the calculation. 
A constant ^3 is included in ^3 to allow the number of 
electrons n; in the /-orbital to be adjusted. 
The total Hamiltonian is given by 



H 



Hk + >^lz K^x - ily) 
) Hk - Xlz 



(2) 



where the first and the second rows (columns) correspond 
to I t) and 11), respectively. The matrix element of l(; 
(C — x, y, z) is specified in quantum mechanics textbooks 
as, for example, {lz)i,m = i{5i,25m,i - 5ui5m,2)- In ac- 
cordance with Refs. [17-19], we set t ^ I (« 0.2 eV), 
t' = 0.1, is = 0.8, 4 = 0.35, and A ~ 0.2. The band 
structure for Sr2Ru04 is reproduced by selecting the 
chemical potential /x and ^3 so that ni = n.2 = "-3 = n/'i. 
The Fermi surfaces for A — 0.2 are shown in the inset 
of Fig. 1. (a,/3)-bands [7-band] are mainly composed of 
{dxz, rfyz)-orbitals [cZajy-orbital] and, due to the symmetry 
of H under k <-> — k, each band is two-fold degenerate 
even if A > 0, in contrast to the Rashba 2DEG model. 
The shape of the Fermi surfaces agrees with that obtained 
in the ARPES measurements, LDA band calculations, 
and dHvA oscillations [20]. 

The 6x6 matrix retarded Green function is given by 
G^( k,w) = [uj+ii-H + iV)-^, where r is the imaginary 
part of the k-independent self-energy due to local impu- 
rities. Since a tiny residual resistivity in Sr2Ru04 sug- 
gests that the impurity potentials are small, we apply the 
Born approximation: F oc — Im^^. (^^(k, 0). If A = 0, 
then F is diagonal; Fj^m oc jiSi^m {h m — 1,2,3). The off- 
diagonal terms will then be negligible if A ^ VFband- 7; 
is the quasiparticle damping rate for the /-orbital, and is 
proportional to the local density of states (LDOS) for the 
/-orbital, pi{0). In Sr2Ru04, 1 < p3(0)//9i(0) - 2.5 since 
the 7-band Fermi surface is very close to the van-Hove 
singular point at (7r,0) [17]. 

The charge current for ^-direction (^ = x, y) is given 
by = — e • dH/dk^ (— e is the electron charge), and 
the CTz-spin current is given by the Hermitian operator: 

independent, they are given by 



(-l/e){J?,s4/2 [4-13]. Since the SOI is k- 
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(3) 



where = dHK / dk^. Also, the /^-orbital current is 
given by jO = {jC,/;}/2 [13]. 

We now calculate the SHC and OHC in the presence 
of local impurities using the linear-response theory. In 
the current model, vq and G are respectively odd and 
even functions with respect to k <-> — k. Therefore, 
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FIG. 1; A-dependence of the SHC and OHC in Sr2Ru04. A 
typical value of A for Ru'''''-ion is ~ 0.2. We set 71 = 72 = 
0.005 and 73/71 = 3. Inset: Fermi surfaces of Sr2Ru04 for 
A = 0.2. 



the current vertex correction (CVC) due to the local im- 
purity potentials in the (extended) Born approximation, 
which is given by AJf oc Y.\.G^JfG^ (M=S,C or O), 
vanishes identically. (In contrast, the CVC is crucial 
in the Rashba 2DEG model since the SOI, which is an 
odd-function with respect to k, gives a k-independent 
anomalous velocity [6].) For this reason, the SHC at 
T = is given by a^y = afy + afj [21], where 
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Here, / and // represent the "Fermi surface term" and 
the "Fermi sea term" . Taking the trace in eq. (4) , (t|^ in 

the present model is straightforwardly given by ax^ 



zl(l) 



zl{2) , 2/(3) , 

(Jxv + cTxv , where 
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where 



dgjdkc^ (C = x,y), Vx = d^i/dkx, Vy = 

d^2/dky, and = det{^i- H" +if). axy"'^ (m = 1, 2, 3) 
is proportional to A™ except for the A-dependence of D^. 
The "Fermi sea term" can also be obtained by inserting 
J^'^ in eq. (3) into eq. (5). We do not give an expression 
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for it since its contribution is very small in the metallic 
state [23] (see Fig. 2 (b)). 

In the same way, the OHC of the Fermi surface term 
O^y and that of the Fermi sea term O^y are respectively 
given by eqs. (4) and (5) by replacing j| with J°. In 
the present model, O^y is given by Oxi}^^ + Oxi}"^^ , where 

= ^ E ^1"- + - ^2)] (9) 

k 
k 

The interorbital velocity = dg/dkQ is termed the 
"anomalous velocity" , which is the origin of the Hall 
effect [23]. Since = — 4i' sinfcj^ cosfc^; has the same 
symmetry as ky, eqs. (6)-(8) are finite even after the k- 
summations. Since t' is much larger than the Rashba SOI 
in semiconductors, a large SHC is realized in the present 
model. Although the sign of t' changes due to the gauge 
transformation \xz) — > — |a^z), both the SHC and OHC 
are invariant if the sign change of {Iq)i,i is taken into 
account [23]. We stress that eqs. (6)-(10) are indepen- 
dent of the rcnormalization factor due to the Coulomb 
interaction: z ~ {1 — dTi{Ld) / duj)~^\i^=Q [23]. 
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FIG. 2: (a) 71-dependence and (b) 73/71 -dependence of the 
|SHC| and |OHCI in Sr2Ru04. \(jly{t' = 0.05)| is larger than 
\(Jxy{'t' ~ 0.1)1, which is understood from the relation |cr|y| oc 
t'A~'^ in the present model, (c) n-dependence of the SHC 
and OHC in the t29-model [Sr2M04]. 

Figure 1 shows the A-dependence of the SHC and OHC 
for the <2£rniodel with Af=Ru (n = 4). Here, we put 



7i = 72 = 0.005 and 73/71 = 3. The Fermi sea terms 
(//) of both the SHC and OHC are neghgible [23] since 
they are about 30 times smaller than the Fermi surface 
terms (/) given in eqs. (6)-(10) (see Fig. 2 (b)). Here, 
1.0 [e/27ra] corresponds to w 670 [he~^ ■ Q^^cm~^] if 
we put the interlayer distance of Sr2A/04, a ~ 6 A. The 
obtained SHC and OHC for typical values of A ~ 0.2 are 
much larger than those in semiconductors [4, 13]. In Fig. 
1, the SHC increases approximately in proportion to A^ 
since a^y^"^^ in eq. (7) gives the dominant contribution 
for 0.05 < A < 0.3. On the other hand, the OHC exists 
even if A = since the orbital current is independent of 
spin. Note that the OHC vanishes if A = in Rashba 
2DEG model since the anomalous velocity arises from 
the k-dependent SOI. 

Figure 2 (a) shows the 71-dependences of the SHC and 
OHC for 73/71 = 3. They are approximately indepen- 
dent of 71 in the low-resistivity region where 71 < 0.1. 
This behavior is typical of the intrinsic Hall effect found 
by Karplus and Luttinger [22]. However, both the SHC 
and OHC start to decrease with 71 in the high-resistivity 
region where 71 > 0.1. This crossover behavior is also 
realized in the intrinsic anomalous Hall effect (AHE) 
[23, 24]: Since the intrinsic Hall effect originates from the 
interband particle-hole excitation induced by the electric 
field [22], the Hall conductivity is proportional to the 
lifetime of the excitation /i/A, where A is the minimum 
band-splitting around the Fermi level [23] . In this model, 
A ~ 0.2. In the high-resistivity region, the SHC de- 
creases drastically since the interband excitation is sup- 
pressed when the quasiparticle lifetime h/j is shorter 
than h/A. Interestingly, the SHC becomes positive for 
71 > 0.3 since axy^^^ (> 0) is dominant over axl^"^^ ■ 

Even in the low-resistivity region, the SHC and OHC 
sensitively depend on the value of 73 /71 as shown in Fig. 
2 (b). \cr^y\ increases significantly with 73/71 since the 
numerator of the second term in eq. (7) and that in 
eq. (8) are proportional to 73. On the other hand, O^y 
is approximately independent of 73/71. Here, we con- 
sider how to increase 73/71 in real systems. According to 
ref. [25], 7-band Fermi surface touches (tt, 0) in La-doped 
Sr2_yLayRu04 for y - 0.23 (n - 4.23). Then, 73/71 > 1 
is expected within the Born approximation. We also con- 
sider Ca-doped Ca2_ySryRu04 {n = 4): According to 
optical measurements [26], a huge mass-enhancement is 
realized in 7-band, whereas a, /?-bands remain light. In 
this situation, 73/71 3> 1 due to inelastic scattering is 
expected at finite temperatures. 

Figure 2 (c) shows the electron filling (n) dependence 
of the SHC and OHC in the t2g-model. For Sr2M04, 
M=Ru, Mo and Rh correspond to n = 4, 2 and 5, re- 
spectively [27, 28]. Here, we correctly account for the 
n-dependence of 73/71 — p3(0)/pi(0). It is thus possible 
to control the values of SHC and OHC by varying the 
electron filling. 
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FIG. 3: "Effective magnetic flux" for j-electron through tri- 
angles of half unit cells in the {d^z, dj;z)-niodel. The effective 
magnetic flux for |-electron is ~(f>o/A. 



We now present an intuitive explanation for the SHE in 
d-electron systems based on the {dxz, dyz)-m.oAc\ [23, 31], 
which is a simplified version of the present model. Here, 
all the matrix elements of Ix and ly vanish identically. 
In Fig. 3, we consider the motion of a j-spin elec- 
tron along a triangle of a half unit cell. An electron 
in the dajz-orbital can transfer to the dj^^-orbital and vice 
versa using the SOI for ].-electron —hXlz, {yz\lz\xz) — 
— {xz\lz\yz) = i. These relations are derived from the 
fact that the (7r/2)-rotation operator about the z-axis is 
given by exp{—i7rlz/2) = cos(7r/2) — i/2sin(7r/2) = —ilz 
for 1^ = 1. By considering the sign of the interorbital 
hopping integral (±i') and the matrix elements of the 
SOI, it can be shown that a clockwise (anti-clockwise) 
motion along any triangle of the half unit cell causes the 
factor -l-i (— i). 

This factor can be interpreted as the Aharonov-Bohm 
phase factor e'^'^^4>/<l>o = hc/e], where = f Adr — 
00 /4 represents the "effective magnetic flux" through the 
triangle of the half unit cell. Since the effective magnetic 
flux for |-electron is opposite in sign, |-electrons and j- 
electrons move in opposite directions. Thus, the effective 
magnetic flux in Fig. 3 gives rise to a SHC of order 
0(A) [eq. (6)], which is analogous to the AHE due to 
spin chirality [29, 30]. (In the i2c,-model, the SHC of 
order O(A^) in eq. (7) is derived from the interorbital 
transition between (xy j) and {xz J.) [or (yz J,)] using 
the SOI twice.) We can also discuss the origin of the 
OHE by considering the motion of an electron with \lz = 
±1) oc ^\xz) +i\yz) in the absence of the SOI. Therefore, 
large SHE and OHE due to such an effective flux will be 
ubiquitous in multiorbital systems. 

We comment that the {dxz,dyz)-mode\ has been ap- 
plied to calculate the large AHE in Ru-oxides [23, 31]. 
Consistent with the theory, paramagnetic compound 
Cai.7Sro.3Ru04 exhibits large AHE under the magnetic 
field [32]. We find that the same mechanism, that is, the 
effective Aharonov-Bohm phase due to c?-orbital degrees 
of freedom, gives rise to the large SHE and AHE. 

In summary, we have studied SHC and OHC in the 



t2g-iiiodel. The SHC obtained is about three times larger 
than the recently observed "giant SHC" in Ft [15], and 
it can be further increased by controlling 73/71. The 
large SHE in our model originates from the "effective 
Aharonov-Bohm phase" induced by the atomic SOI and 
the interorbital hopping integral. Moreover, huge OHC 
in this model will enable us to control the atomic d-orbital 
state by applying an electric filed. The present calcula- 
tion predicts that "giant SHE and OHE" will be ubiq- 
uitous in d-electron systems. Using the theoretical tech- 
nique developed here, SHC in various transition metals 
such as Ft had recently been investigated [33]. 
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